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We investigate a string-inspired dark energy scenario featuring a scalar field with a coupling to 
the Gauss-Bonnet invariant. Such coupling can trigger the onset of late dark energy domination 
after a scaling matter era. The universe may then cross the phantom divide and perhaps also exit 
from the acceleration. We discuss extensively the cosmological and astrophysical implications of 
the coupled scalar field. Data from the Solar system, supernovae la, cosmic microwave background 
radiation, large scale structure and big bang nucleosynthesis is used to constrain the parameters of 
the model. A good Newtonian limit may require to fix the coupling. With all the data combined, 
there appears to be some tension with the nucleosynthesis bound, and the baryon oscillation scale 
seems to strongly disfavor the model. These possible problems might be overcome in more elaborate 
models. In addition, the validity of these constraints in the present context is not strictly established. 
Evolution of fiuctuations in the scalar field and their impact to clustering of matter is studied in detail 
and more model-independently. Small scale limit is derived for the perturbations and their stability 
is addressed. A divergence is found and discussed. The general equations for scalar perturbations are 
also presented and solved numerically, confirming that the Gauss-Bonnet coupling can be compatible 
with the observed spectrum of cosmic microwave background radiation as well as the matter power 
spectrum inferred from large scale surveys. 
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I. INTRODUCTION 

General relativity predicts singularities. Therefore, in 
spite of its being a highly successful as a classical theory 
of gravitation, its modification seems inevitable at high 
energy scales. The low energy action could then also fea- 
ture corrections to the Einstein-Hilbert term from addi- 
tional fields and curvature invariants. Among the possi- 
ble corrections, a particular combination of the quadratic 
Riemann invariants, the Gauss-Bonnet term, is of spe- 
cial interest. It appears in extensions of gravitational 
actions, whether motivated by the form of most general 
of a scalar-tensor theory, uniqueness of the Lagrangian in 
higher dimensions or the leading order corrections from 
string theory. An interesting possibility then arises that 
the modifications could change the way the universe itself 
gravitates at large scales. 

Indeed, the current cosmological observations cannot 
be explained by the standard model of particle physics 
and general relativity. Precise cosmological experiments 
have confirmed the standard Big Bang scenario of a flat 
universe undergoing an inflation in its earliest stages, 
where the perturbations are generated that eventually 
form into structure in matter. Most of this matter must 
be non-baryonic, dark matter. Even more curiously, the 
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universe has presently entered into another period of ac- 
celeration. Such a result is inferred from observations of 
extra-galactic supernovae of type la (SNela) and is 

independently supported by the cosmic microwave back- 
ground radiation (CMBR) [3| and large scale structure Q 
data. It seems that some dark energy, with its negative 
pressure that speeds up the universal expansion, domi- 
nates the density of the universe Q. This concordance 
model agrees very well with the astrophysical data, but 
features inflation, dark matter and dark energy as phe- 
nomenological ingredients of undisclosed nature. Perhaps 
dark energy is the most enigmatic of the latter two, as 
at least there are reports of tentative discoveries of dark 
matter Q. 

The problem of dark energy reduces to the questions 
1) what is the culprit (perhaps in fundamental physics) 
for the acceleration, and 2) why did the speed-up begin 
just at the present stage of the cosmological evolution. In 
recent years dark energy and its cosmological and astro- 
physical signatures has been addressed in many papers 
considering both modifications of the energy-momentum 
tensor, with the inclusion of scalar fields [3, H, Qii] or 
imperfect fluids 11 , |lj and of gravitational physics at 
large scales [H, El 111, 13 • 

The latter approach has 
proven to be perhaps surprisingly difficult. Extending 
the action with curvature invariants rather generically 
results in fourth-order gravity, which is problematical 
because of the implied instabilities An f(R) exten- 
sion, although not generically, can be stablepSj. How- 
ever, though some f{R) -type modifications might be 
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compatible with the acceleration following a standard 
matter dominated phase [l^ the simplest modifica- 
tions as proposed this far seem not to be able to gener- 
ate a viable cosmology [2lj. On the other hand, within 
the Palatini formulation, the corresponding extensions 
result in second order field equations, but it turns out 
that the observationally allowed modifications are then 
practicall y in dis ting uishable from a cosmological con- 
stant [nin, m if . 

Scalar-tensor theories of gravity are interesting alter- 
natives to the concordance model and seem to have po- 
tential to provide a linka ge b etween the acceleration and 
fundamental physics [26l l27l l28l [29| . They have a desir- 
able feature which is their quasi-linearity: the property 
that the highest derivatives of the metric appear in the 
field equations only linearly, so as to make the theory 
ghost free. Interestingly, there is a particular combina- 
tion of the curvature squared terms with such behavior, 
known as the Gauss-Bonnet (GB) integrand. It is con- 
structed from the metric as 
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As mentioned in the beginning, this term appears fre- 
quently in attempts at a quantum gravity, especially in 
stringy set-ups. In fact, all versions of string theory (ex- 
cept Type II) in 10 dimensions {D = 10) include this 
term as the leading order a' correction [sO, IM. InD = 4: 
the Gauss-Bonnet term is a topological invariant. It's ap- 
pearance alone in the action can then be neglected as a 
total divergence (though even then this topological term 
can have interesting role in view of the boundary terms 
and regularization of the D = 4 action in asymptotically 
AdS spacetimesfH, H^). However, if coupled, its pres- 
ence may lead to contributions to the field equations also 
in _D = 4. The low energy string action typically features 
scalar fields with such couplings. Namely the moduli, 
associated to the internal geometry of the hidden dimen- 
sions, become non-minimally coupled to the curvature 
terms. Hence, in a general _D = 4 scalar-tensor theory 
one might couple the scalar field non-minimally not only 
to the Ricci curvature R but to the Gauss-Bonnet in- 
variant Rqs as well, and in Z? > 4 it is necessary to 
include this term in order to preserve the uniqueness of 
the gravitational action. 

Such a coupling could thus be useful in modeling both 
the early inflation and the late acceleration [sj, [s^, [H, 
\3T\ . [ssj . Here we will concentrate on the post-inflationary 
epochs in such a universe, and more specifically, within 
the scenario we recently studied '39*1 . There the scalar 
field lives in an exponential potential, and the coupling 
can likewise be exponential. This kind of model was orig- 
inally proposed by Nojiri, Odintsov and Sasaki^^l- More 
recently it was shown that if the coupling grows steeper 
than the potential decays with the field, acceleration can 
occur with a canonic scalar, after a transition from a 
scaling erajs^. Then the interaction with the Gauss- 
Bonnet curvature term causes the universe to enter from 
a scaling matter era to an accelerating era. Existence 



of scaling solutions with more general parameterizations 
has then been taken under investigation 4l| , noting that 
requiring exact scaling is otherwise possible only with a 
tachyon field which is non-minimally coupled with mat- 
ter. Other dark energy solutions in a broad variety of 
second order string cosmologies, taking into account both 
coupled and uncoupled matter, have been explored also 
previously [42, . .43] , higher order terms have been also 
incorporated and a method to reconstruct the coupling 
and the potential has been found [i^. If the scalar field is 
non-dynamical, the model is equivalent with a modified 
Gauss-Bonnet gravity, featuring a function of Rq^ the 
action [ii,!!!]. 

The coupling mechanism can also momentarily push 
the universe to a phantom era, until the scalar field po- 
tential begins to dominate. Therefore it is possible to 
get a very good accordance with the SNela data. Dark 
energy models other than of the Gauss-Bonnet form but 
featuring transition to phantom expansion and having 
possible origin in string theory have been considered also 
147', 'isj . In addition, one should clearly distinguish the 
models investigated here from those featuring a coupling 
to the Ricci scalar [4^, since for them the so-called R- 
boost [5C| plays a role in the early universe, while the po- 
tential has to be shallow in order to have accelerating be- 
haviour at the present. On the other hand, one can cou- 
ple the self-tuning scaling field to matter [5l|, [H, , and 
thereby indeed set the acceleration ongoingand thereby 
perhaps alleviate the coincidence problem |M . [ssj . How- 
ever, implications for the large-scale structure formation 
seem to make this possibility problematical [H, [s^l • Here 
we will not find the same, since modifications to the 
growth of matter perturbations are much less drastic as 
the scalar field remains minimally coupled to the matter 
sector. 

In this paper, besides investigating the cosmological 
evolution of such models, we also derive quantitative 
constraints on the scenario of Gauss-Bonnet dark energy 
mentioned above. Our aim is to present in more detail 
this scenario and to generally uncover in more depth what 
are the possible effects of the Gauss-Bonnet combination 
on dark energy cosmologies, in order to assess whether 
they could be compatible with present days observations, 
and to estimate the amount of fine-tuning that this might 
require. To this end the the model in Ref.[39| is sub- 
jected to further constraints. In addition to tightening 
the limits for this particular parameterization, we will 
discuss in more detail and generality the phenomenology 
of Gauss-Bonnet dark energy at cosmological, astrophys- 
ical and Solar system scales during the cosmological evo- 
lution. Especially, we consider the impact of the coupling 
to the CMBR and large scale structure. In view of the 
remarks above, this is a crucial aspect of a dark energy 
model with any non-minimal couplings or modifications 
of gravity. Since the class of models studied here features 
both and these can lead to completely different predic- 
tions for the expansion rate of the Universe than without 
the Rqq term, it is rather non-trivial to find that the 
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perturbation evolution in these models is in some sense 
only modestly altered. Namely, the shape of the mat- 
ter power spectrum is retained during the presence of 
the Gauss-Bonnet gravity and the effects one expects to 
see in the large multipoles of the CMBR spectrum are 
typically small. Still, there are distinguishable features 
peculiar to these models that can lead to essential con- 
straints for the model parameters. 

The article is organized as follows. The cosmological 
dynamics for these models is described in Section |TT] for 
both the background expansion and for linear pertur- 
bations. In Section Uni we discuss the implications and 
derive constraints of the Gauss-Bonnet coupling to nu- 
cleosynthesis, Solar system, large scale structure, CMBR 
and SNela. Section IIVI contains our conclusions. Some 
technical details are confined to the Appendixes. 



II. COSMOLOGICAL DYNAMICS 



Before writing down the cosmological equations, let us 
briefly discuss the model and its motivation. A peculiar 
property of the string effective action is the presence of 
scalar fields and couplings which are field dependent, and 
thus in principle space-time dependent. The scalar fields 
are moduli associated to geometrical properties of com- 
pactified extra dimensions. Thus one could consider mul- 
tiple scalar fields, representing the dilaton and various 
other moduli, but here we stick to one field </>. The rea- 
son is that the dilaton couples to the Ricci curvature, and 
this would lead to variations of the Newton's constant. 
When going from the string frame to the Einstein frame, 
the coupling to R is transformed into a non-minimal in- 
teraction with matter. This in turn would lead to viola- 
tions of the equivalence principle. Thus, in general, also 
the matter Lagrangian L„i could be non-minimally cou- 
pled to the scalar field. There are tight constraints from 
observations to such effects [5^ (see however, [sol f60|). 
However, it is (usually) only the dilaton field which in 
heterotic string theory acquires the geometrical coupling 
to R and thus enters in the conformal factor and the 
matter Lagrangian in the Einstein frame [6l[. If the con- 
formal factor is not (nearly) trivial, it can evolve so, for 
instance due to the so called least coupling principle 6z| . 
In the present article we will make the usual assump- 
tion in cosmology, that possible non-minimal couplings 
to the matter sector are negligible, since our specific pur- 
pose here is to study the novel features resulting from the 
Gauss-Bonnet curvature interaction. Thus we identify 
as a run-away modulus without direct matter couplings. 
Then the equivalence principle is automatically satisfied, 
though gravitational dynamics are modified due to the 
modulus-dependent loop corrections. 

Thus we can write the action for the system to consider 
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here n = {8t:G)~^^^ , G being the Newton's constant. 
The 7 could be a function of 4> as well, and there could be 
additional kinetic terms. In the following we will mostly 
consider canonic scalar field and thus set 7=1 (which 
can be achieved by a redefinition of the field for any con- 
stant 7 > 0). V{(j)) is the field potential which could 
result from nonperturbative effects. In four dimensions 
the GB term makes no contribution if f{4>) = const. 
However, it is natural to consider /(0) to be dynami- 
cal. This follows, for example, from the one-loop cor- 
rected string effective action [61I, [6^ , where the function 
/((/)) = (T — (5^((/)): the coupling a may be related to string 
coupling gs via a ~ l/<?s j ^'^^ numerical coefficient S 
typically depends on the massless spectrum of every par- 
ticular model [131 • It thus seems that the dilaton might 
couple already at the string tree-level next to leading or- 
der expansion in the inverse string tension a', and the 
modulus functions ^{4>) are non-trivial at the one-loop 
order. 

In the numerical examples we will adopt an exponential 
form for the potential and the coupling, 



y(^) = Vbe-^'^'^/^, / = /oe"'^'^/^. 



(3) 



On one hand, the nonperturbative effects from instan- 
tons or gaugino condensation typically result in an ex- 
ponential potential. On the other hand, there is also 
phenomenological motivation for this, since, besides its 
being a simple choice, an exponential potentials allows to 
consider scaling solutions in cosmology. An exponential 
field-dependence for coupling is a reasonable assumption 
in supergravity actions. For massless dilaton one in fact 
has f{(f>) — Y^Cne'-^~^^'^ [11]. For another instance, a 
known example from heterotic string theory [63 | yields 
for the modulus coupling, in terms of the Dedekind func- 
tion ry, the form /(0) - log [2e"'='^r?4(ie''"'*)], which be- 
haves like Eq.([3]) to a good approximation. Again, it is 
a simple choice for /, introducing but one extra param- 
eter a. Such minimalism is practical when considering 
a phenomenological model of dark energy, though from 
a particle physics point of view one might expect that 
corrections to the Eqs.Q come into play, in particular 
when the Gauss-Bonnet term begins participate in cos- 
mological dynamics in the late time in our model. 

In the presence of the coupling /, the field equations 
can be written as 
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where the two first terms in the right hand side are the 
energy-momentum tensors for matter and the scalar field, 
respectively. The curvature corrections resulting from 
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taking into account the Gauss-Bonnet contribution in- 
volve only terms proportional to derivatives of A 
partial integration of the coupled term in action ([2]) gives 
a vanishing boundary term but in addition, because of 
the field dependent interaction, a contribution involving 
the integral of Rq^ which is of the first order by con- 
struction. Explicitly, one obtains 
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The equation of motion for the scalar field reads 

^□^_r(</>)-/'(0)i?^s=O. (6) 

Thus the field lives in an effective potential given by 
V{4)) + f{4>)R%B- Matter, since minimally coupled, is 
conserved as usually [11], = 0. The general covari- 

ant expression for curvature corrections Eq. ([5]) is a rather 
complicated combination of the derivatives of the cou- 
pling and their contractions with the Riemann and Ricci 
objects, but due to the geometric properties of these cor- 
rections, they vanish in constant curvature spacetimes, 
and also otherwise can assume rather tractable forms, as 
we will see in the following. 

A. Background 

We consider a flat, homogeneous and isotropic back- 
ground universe with scale factor a{t) in the Friedmann- 
Robertson- Walker (FRW) metric 

ds^ = -dt^ + a^{t)6^jdx'dx^ . (7) 

The action ^ then yields the Friedmann equation as the 
t-t component of Eq.(|4]), 



where an overdot denotes derivative with respect to the 
cosmic time t, H = d/a is the Hubble rate and Pm repre- 
sents the matter component. The Klein-Gordon equation 
Eq.® reads 
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7(0 + 3H<i)) + y + n<j,)RlB = 0, 
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where the Gauss-Bonnet invariant is R%g — 2AH^ [H + 
H^). It will be convenient to work with the dimensionless 
variables defined as 
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In addition it is useful to define, in analogy with fim, 
the relative contributions from the scalar field and the 
Gauss-Bonnet correction as 



JlP ~ 3 ' 



f 



1 ^m. 



^l. (11) 



Then 



Weff = W„iQ.„i + + WfQ.f = - 1 (12) 

is the total effective equation of state in the sense that 
the Universe expands as if dominated by a fiuid with 
this relation between its pressure and energy. One notes 
that the Gauss-Bonnet term can be written as R^g = 
— 12H*{3'Wef f + 1), and is thus negative if and only if the 
scale factor is decelerating. 

With the aim to understand the behaviour of the back- 
ground cosmology for the model we perform a dynamical 
system analysis. In terms of the dimensionless variables 
(jTO|) . the complete dynamical system is then given by [s^ 
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-3 + 72;^ + y + 3/i + 3ilm, 
2e + 3 + jx'^ — y — ^' — (e + 2)^ + Swm^m, 
27 [xx' + x^{e + 3)] + y' + 2ye + 3^i{e + 1), 
2xx' + y' + 3Ai' + 3f7'„, 

-3(1 + W„i)flm 



2en,, 



(x' jx -I- 2e ax)pL, 
-{Xx + 2e)y, 



(13) 



where prime means a derivative with respect to log(a) 
and Wm is the equation of state of the background. The 
first five equations hold generally for an action like 
whereas the two last equations encode the information 
about the specific model ([3]). For simplicity, we set from 
now on 7 = 1, so the scalar field kinetic term is canonical. 

We find several fixed points for the system (|13p. char- 
acterized by x' ~ y' = fi' = 0: 

• A: {x,y,fi) = (0,0,0). This is an unstable point 
which would correspond to domination of the back- 
ground fluid with Weff = Wm and completely van- 
ishing contribution of the scalar field. 

• B: {x,y,ii) = (±V3, 0, 0). In this case, the ki- 
netic energy of the scalar field dominates. Then 
Weff = 1. The kination phase is not stable, since 
the kinetic energy will always redshift away faster 
than the other contributions to the energy density. 
This solution with > is a saddle point (it has 
one positive eigenvalue, implying that the solution 
attracts from some direction while repelling from 
some other), whenever a < 2^/3 or A, a > 2^3. 

• C: (x, ?;,/i) = (0,3,0). In this fixed point, the po- 
tential of the scalar field dominates. Thus this is a 
de Sitter solution with Weff = — 1. In the absence 
of the coupling, all other contributions to the en- 
ergy density have larger effective equations of state, 
and the solution is stable. With the coupling, sta- 
bility condition is simply a > A. 

• D: {x,y,p,) = (A/2,3 - AV4,0). This fixed point 
corresponds to scalar field domination with We// — 
A^/6 — 1. It does not exist when A > 2\/3 and is 
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thus irrelevant to us here. When a = A, the Gauss- 
Bonnet term can reduce w^ff of this solution(4l|. 

• E: ix,y,fi) = (f(H- «;™),^(l - zz;^),0). This 
is the well-known scaling solution, where a scalar 
field with exponential potential mimics exactly the 
background equation of state Wm- This fixed point 
is a stable spiral when A > a, \/6(l + Wm)- It is, 
however, a saddle point when A, a < \/6(l + w) or 
a > A, a/6(1 + Wm). Also this solution exists in a 
modified form in the special case a = A [4l| . 

.F: = (f(l + ^„0,0,.^ ^^:T+\^y"^^ ). 

This is an a new scaling solution, in which field 
has run to large values and the the potential can 
be neglected. Then 57^ is nonzero and the effective 
equation of state in this case is again just Wm- The 
condition for the existence of this fixed point is that 
either Wm > 7/3 or = 3(1-1- Wm)^(7—3wm)/(l + 
3wm)- Then it might be a saddle point, but it is 
generically unstable. Thus this possibility of the ki- 
netic energy together with the Gauss-Bonnet con- 
tribution scaling like matter seems not to be use- 
ful for cosmological applications, as it will not be 
reached by dynamical means. Even if one sets the 
field to this solution as an initial condition during 
radiation domination, the Gauss-Bonnet contribu- 
tion will drop away when transition to matter dom- 
ination occurs. 

• G: There exists also a solution where only the 
Gauss-Bonnet term and the kinetic term of the 
scalar field survive while y — flm = (see Ap- 
pendix This solution is complicated and of no 
cosmological interest to us here. 

From this we can see that the standard tracking be- 
haviour of exponential quintessence is available whenever 
the coupling term is negligible. This tracker solution has 
been shown to exist for extremely wide range of initial 
conditions. Unfortunately, while in the tracking regime, 
the scalar field equation of state equals exactly the back- 
ground w, and thus this solution cannot account for the 
accelerating universe. 

If, however, the coupling becomes significant at late 
times, the situation will change. Since for the fixed point 
E we have ~ p ~ ^-3(1+™^)^ ^^le last term in the 
Friedmann equation ([8]) scales like pf = H^f'{(l))(j) ^ 
Q-^3(i+»„)(2-a/A)_ rpj^jg follows from the tracking be- 
haviour of the scalar field; since (p' = 3(1 -I- w,,n)/X, we 
have that (/> = 0o + 3(1 -I- log(a)/A, which implies 
f{(j)) - a^o:{i+w^)/\_ rpj^^g gj^^ ^^^^ effective en- 
ergy density due to the presence Gauss-Bonnet term, p f , 
dilutes slower than the energy density due to ordinary 
matter, p™, if and only if a > A. We have confirmed this 
simple result by numerically integrating our system (jl3p 
(see Appendix VK\ for details). Then we found that as the 
coupling begins to affect the evolution of the field, it will 
always be passed to the fixed point C from the saddle 
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FIG. 1: A phase portrait of the model. The saddle point 
E attracts the field to the scaling regime but later always 
passes it, along the same track, to the fixed point C which 
corresponds to an accelerated expansion. The requirement for 
the fixed point C to be reached and be stable is a > A. For 
the acceleration to begin only after a sufficiently long scaling 
matter era, the scale /o has to bet set suitably. 



point E as depicted in FIG[T] Hence the universe is ap- 
proaching a de Sitter phase and an acceleration occurs 
as observations indicate. 

Furthermore, during the transient epoch between mat- 
ter and scalar field domination, the universe can enter 
into a transient phantom stage (with Weff < 1), which is 
possible due to the pf term in the Friedmann equation. 
Even if Weff > —1, a result Weff < —^4, would imply 
effective phantom dark energy when is interpreted as 
due to an uncoupled dark energy component. We plot 
two typical examples of such evolution in FIG. [5] with 
A = 4 or A = 8 and a = 20. Note that Q„i + > 1 is 
possible, since the (effective) Gauss-Bonnet energy den- 
sity pf can be negative when the field momentarily rolls 
backwards. This can occur since, as noted previously, 
when the universe begins to accelerate the Gauss-Bonnet 
term flips its sign, and might overturn the slope of the 
effective potential. 

Having now found the background expansion in the 
model, we can check what parameter values are needed 
to produce it. We sec from Fig.(l2|) that all the variables 
X, y and p are roughly of the order of one today. Then 
we have that 

if we set /o of order one as expected in string theory and 
recall that {k^H)'^ « lO^^^. Note that though f{(f>o) has 
to be very large, because of the exponential amplification 
oi f{4>) we can here have avoid introducing huge numbers 
into the Lagrangian ([2]). Then it is seen that the present 
value of the scalar field is about k4)q ~ 390/a. It seems 
that typically our models the field has run to perhaps a 
couple of dozen Planck masses. Then the scale of the 
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FIG. 2: The fractional energy densities for matter, flrn (dash- 
dotted line), the scalar field fl^ (dashed line) and Gauss- 
Bonnet term, Qf (dotted line). The solid line is the total 
equation of state Wef / . The upper panel is for A = 4 and the 
lower panel for A = 8. For both cases — 0.35 and a — 20. 
The transient phantom era in the upper plot is caused by the 
dynamics of the coupling. 



potential can be deduced from 



1. 



(15) 



In the case X — a the potential scale Vq would be of 
the order of Mp, where Mp — 1/k is the Planck mass. 
The mass of the field at present, = \/—VQ{(f>o), on 
the other hand, will turn out to be very small, as usual 
with quintessence fields. The potential energy is expo- 
nentially suppressed, while the coupling is exponentially 
enhanced as a function of the field; this allows to con- 
sider natural magnitudes for both simultaneously. How- 
ever, for the same reason these estimates are also very 
sensitively dependent on the particular values of the pa- 
rameters. Since in the models we consider here a > A, 
the potential scale is in fact usually much less than Mp. 
According to the very rough estimates and 
if /o ~ 1, then a could be allowed to be only slightly 
less than A, for the potential scale to be of the order of, 
say Vf^^^ - (to3/2Mp)^/2 - IQiO-lO" GeV, which cor- 
responds to the effective potential for low energy super- 
symmetry breaking, with the soft supersymmetry mass 



scale of about m3/2 = 10^ GeV. Nevertheless, this sug- 
gests that the scales of the model could be obtained from 
more fundamental physics. 



B. Linear Perturbations 

The line-element Eq.([7|) generalizes in the perturbed 
FRW spacetime to 



- (1 + 2ipA) dt^ ~ 2a{t)l3^tdtdx'' 
a^[t) [5ij (1 + 2ip) + ipi,j] dx'-dx' . 



(16) 



We consider variables in the Fourier space. The trans- 
formation is simple since at linear order each fc-mode 
evolves independently. We do not consider vector per- 
turbations since their evolution is unmodified(67| . and we 
will only briefly comment on tensor perturbations. We 
characterize the scalar perturbations in a general gauge 
by the four variables (pA, P, <f, ip- Some of these degrees of 
freedom are due to arbitrariness in separating the back- 
ground from the perturbations. One can deal with these 
gauge degrees of freedom by noting that the homogene- 
ity and isotropy of the background space implies invari- 
ance of all physical quantities under purely spatial gauge 
transformations [68]. Therefore one can trade P and ip to 
the shear perturbation 



(17) 



Since both [3 and ip vary under spatial gauge transfor- 
mation, they appear only through the spatially invariant 
linear combination x i^i relevant equations. In addi- 
tion, one can define the perturbed expansion scalar 



(18) 



Using these variables in the so called gauge-ready formal- 
ism is useful in studying generalized gravity [22, [63, |63]. 
The general equations for scalar perturbations are given 
in Appendix IB] From here on we set (87rG')~^/^ to unity 
to avoid confusion since it was also denoted by n. 



1. Small scales 

Because the Gauss-Bonnet interaction becomes here 
dynamically important at late times, let us for now ne- 
glect the contribution from radiation and consider the 
matter dominated and subsequent epochs, when Wm — 
= Tim = 0. Dropping the subscripts w, the continuity 
equations (jB8p and (|B9|) then read 

^-Hv + -(pA. (19) 
a 



— V 

a 



K - 3H(pA, 



We choose to work in a synchronous gauge, defined by 
ipA = 0. One can then see from Eg. ([TO]) that it is con- 
venient and legitimate to set u = to fix the remain- 
ing gauge mode. Then Eq. (fTO|) tells that k = 6. Thus 
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already three variables are eliminated from the system. 
The momentum constraint (|B3p then allows to relate the 
evolution of metric potential ip to the fluctuations in the 
scalar field by 



+ H^f'6(j). (20) 



Note that in the conventional synchronous gauge nota- 
tion [t^I k = —h/2 and 77 = — cp, from which follows that 
the metric shear x — ^a^{h/2 + 3ri)/k'^. The evolution 
of this X is governed by 



(1 - 8Hf)x +[H- 8iHf + Hf + H^f)]x 



(21) 



Now we turn to consider small scales in order to find 
approximations for the perturbation evolution. Then 
carefully neglecting the terms that are subdominant at 
the large-fc limit yields the energy constraint (|B2p in the 
form 



(1 - 8ff/)^ -{H- 12H^f )x - iH^SfScf>.{22) 



Equations ^U\\ and (|2ip imply that as usually, here the 
metric fluctuations ip and x (oi' v) ^re small at large k, 
since their only sources are the fluctuations of the scalar 
field. However, from Eg. p^ we do not expect this to 
be the case for k (or h). The large part of the spatially 
gauge-invariant variable k is separated, by the definition 
p8)) . into the gradient of the metric shear Xi which we 
have kept in Eq. . We don't either drop the scalar field 
fluctuations or their derivatives, since their gradients can 
influence the matter inhomogeneities. Since the field is 
very light, <^ fc^ is easily satisfied at the scales we 
consider, and the potential term could be dropped from 
the perturbed Klein-Gordon equation. Finally, we have 
the small-scale limit of the Raychaudhuri equation (jB5[) . 



a?p5 



(1 - 8ff/)x + m'f - AH^f + 8Hf)x 



8{f~Hf)p^4{H^ + 2H)f6 



(23) 



With some algebra, it is now possible to deduce the evo- 
lution equation for the matter overdensity d. One arrives 
afi 



2HS = AirCpS. 



(24) 



The effective gravitational constant G* seen by the mat- 
ter inhomogeneities depends rather non-trivially on the 



^ Assuming the BirkhofT theorem (or Jebsen-BirkhofT theorem (23ll ) 
and deriving the evolution equation for spherical overdensities 
at subhorizon scales as suggested Ref . l7l| for some different 
modified gravity models than the present, would give G«~^ = 
2[2e(l -I- e) -I- e']/(3f2m). The reason for discrepancy is that the 
Jebsen-Birkhoff theorem is not respected in the model. 



evolution of the background quantities, and in terms of 
our dimensionless variables defined in Eq.ljTU]) it can be 
written as 



G* =4 



-x^ + m'(1 + e)' + x^[2{l + £)(/i -l) + y] 
x2[4 + p{5fi - 8)] - m2[6(1 + e)ip-l)+ y] 



(25) 



We plot G^,/G in the upper panel of Fig. [3] for selection 
of parameter values. It is clear that at high redshifts, 
when the scalar field and so coupling term is negligible, 
all the models reduce to the standard value of G. 

When the coupling can be neglected (/',/" ~ 0), G* 
reduces to unity. This is understandable since standard 
uncoupled quintessence does not cluster at small scales, 
and the evolution equation for linear growth stays un- 
modified. On the other hand, in the slow-roll limit where 
(j) and (f) can be neglected we have 



l + 32f2(-2g2 + p)2 
l + 32i72/'2(3i72_2p)' 



(26) 



In the case that the coupling is subdominant, so that 
H'^f',H^f" - e, one has 



V'H) 



(27) 



G, = 1-8 

+ 32f'[^^r~44»PfH 

- ^^f'H + f'{H^-pf + 2^^f'{2H^ + p) , 



where the first and second square brackets respectively 
embrace corrections of the order 0{e) and 0{e^). 

In any case, we can here make the important conclu- 
sion from Ea. ([M)) that the shape of the matter power 
spectrum is the same as for ACDM cosmology. The su- 
perhorizon scales, where our approximation breaks down 
and the density perturbation becomes gauge dependent, 
are not effeciently probed by the present large scale struc- 
ture surveys. As we have shown that at subhorizon scales 
the growth rate of structure is the same at all scales, 
which is the also case when the acceleration is driven by 
vacuum energy, we cannot distinguish between these two 
very different scenarios by comparing the shape of matter 
power spectrum (assuming of course that the primordial 
spectrum is the same in both of the cases) . Therefore the 
primary constraints arising from the matter power spec- 
trum could be deduced from the overall normalization. 
Note that these conclusions are model independent, the 
only assumption being that the scalar field is not very 
massive. We did not make specific assumptions on the 
coupling or the potential or assume that flm dominates 
in order to deduce these results. The constraints from 
large scale structure are still of course relevant, since it 
is possible that for some models the (scale-independent) 
growth rate is heavily modified and the normalization of 
the power spectrum can be used to exclude such cases. 
On the other hand, the perturbations at large scales will 
affect the integrated Sachs- Wolfe effect of the CMBR, 
which might constrain models although the cosmic vari- 
ance weakens the significance of the low CMBR multi- 
poles. 
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FIG. 3: Evolution of effective gravitational constant G* (top 
panel) and of the dimensionless growth rate dlog Sm / {dlog a) 
(lower panel) as functions of redshift. We see that a diver- 
gence is possible for some parameter combinations. The thick 
lines are numerical solutions of the full linearized equations, 
and the thin dotted lines in the lower panel are solutions to 
the approximative equation (|24p . The agreement is excellent 
in most of the parameter space, though in the two extreme 
cases depicted here (dashed line corresponding to very slow 
transition, the dash-dotted line corresponding to instability) 
deviation is visible. 



2. Numerical solutions 

To investigate the cosmological effects of the Gauss- 
Bonnet coupling in more detail, we have numerically in- 
tegrated the full perturbations equations^ and computed 
the full matter power and CMBR spectra for the exam- 
ple model presented in the previous subsection. Then we 
have to fix several other parameters besides the il'^, A 
and a. Since these parameters govern the background 
expansion and thus fix the distance to the last scatter- 
ing surface, they determine the locations of the peaks 



^ In numerical calculations we use a modified version of the public 
CAMB code[7^, see |http: / / camlT. info / Our numerical imple- 



in the CMBR spectrum, as will be seen in next sec- 
tion. However, the relative peak heights depend on the 
amount of baryonic matter. To get results that are 
(roughly) compatible with the WMAP observations we 
fix rijj/ri^ — 1.716 for all the cases shown in this section. 
For simplicity, we set the scalar spectral index ns to one 
and the optical depth r to zero and keep the present Hub- 
ble constant h fixed to 0.7. The normalization we use 
for each model is such that the first peak in the CMBR 
spectrum matches with WMAP observations. For com- 
parison, we plot also the concordance ACDM model with 
its parameters set according to the best-fit obtained by 
combining the WMAP and SDSS data [S]. 

Firstly we can confirm the considerations based on 
appropriate approximations and to compute the linear 
growth rate found in the previous subsection. The di- 
mensionless quantity F — (log£(5))' can be introduced to 
characterize the growth rate[23j. We have checked that 
its value is now independent of scale in the matter dom- 
inated era, when k > O.Q2h^^. Note that there 73] 



Fmd — 1 



(28) 



mentation is discussed in the Append: 



ixes. 



which is the small scale solution to Eq. ([24j) when the 
Gauss-Bonnet terms can be neglected. The subsequent 
evolution of this quantity is plotted for various parameter 
choices in the lower pannel of FIG. [3l showing that this 
evolution is indeed reproduced to high accuracy by the 
simple second order differential equation Eq. ((24|) . Dis- 
crepancy with the numerical solutions to the full lin- 
earized equations appear only for some extreme cases. 

The effects on the CMBR and matter power spectra 
of varying the model parameters is shown in FIGS. 

Note that though we include the error bars in the 
figure, they only roughly indicate how the observations 
constrain the models, since likelihoods should be com- 
puted using the window functions that depend on each 
model. The CMBR error bars are from the three year 
WMAP data 3 and the error bars for the matter power 
spectra are provided by SDSS [2| . A complete likelihood 
analysis taking into account all these data would require 
exploring a vast parameter space, which is beyond the 
scope of the present study. Here we are rather interested 
in study how the model qualitatively differs from the con- 
cordance cosmology when the inhomogeneous evolution 
is considered. In the next Section we will see that in fact 
it is enough to impose tight constraints on the model 
parameters considering only the background expansion. 

FIG. ^ seems to indicate that the model favours mat- 
ter densities about ~ 0.4. It seems clear that we 
cannot do without dark matter within this framework. 
Low matter densities result in large relative contribution 
from the Gauss-Bonnet term (this means very negative 
effective equation of state, as shown before), and leads 
to, in addition to a modification of the peak structure in 
CMB, a large ISW effect. Requiring viable normalization 
for both the CMBR and the matter power spectrum can 



9 




0.005 0.01 



0.05 0.1 
k(h Mpc^') 




0.005 0.01 



0.05 0.1 
k(h Mpc"') 



FIG. 4: The effect of matter density on the CMBR and matter 
power spectra. Here A = 6 and a = 20. Dotted hnes are for 
SlJii = 0.3, dashed hne for SI^ti = 0.4), and dash-dotted for 
n° = 0.5. The soUd Une is ACDM model. 



FIG. 5: The effect of the potential slope on the CMBR and 
matter power spectra. Here Q,^ = 0.4 and a = 20. Dotted 
lines are for A = 4.5, dashed line for A = 6.0), and dash-dotted 
for A = 8.0. The solid line is ACDM model. 



3. Stability of perturbations and a possible end of 
acceleration 



also significantly constrain the allowed matter density. In 
FIGs. (O and (O we show how the cosmological predic- 
tions are changed when the slope of the potential or of the 
coupling are varied. The main imprint from different po- 
tential slopes A seems to be in the normalization. For low 
values of A, there is significant contribution of the scalar 
field during the scaling matter era. This slows down the 
rate of growth of matter inhomogeneities, see Eg. ipS)) . 
Hence the fact that there is less structure nowadays than 
for larger A is not a consequence of the Gauss-Bonnet 
modification, but rather an effect of the presence of the 
scalar field in the earlier scaling era. Finally, in FIG. 
^ we see that the strenght of the coupling a might be 
more difficult to deduce from these data. With steep cou- 
pling slopes, the scalar field domination takes place more 
rapidly and with more negative w^ff, which can some- 
what amplify the ISW effect. The contrary happens for 
smaller a. 



The effective gravitational constant Eq. (|^5)) can in 
principle diverge . Indeed, we have found that in the 
model studied here G* typically diverges in the future, 
see FIG. [31 For low matter densities or large a this can 
happen even before a = 1, as happens in an example 
plotted in FIG. ([3]). Such a case would clearly be ruled 
out. However, in the model parameter combinations that 
would lead to this divergence of perturbations before the 
present day, though may fit the SNela or other individ- 
ual data, are not usually compatible with the combined 
constraints (these will be presented in the next section). 

Matter perturbations will at least for a while grow ex- 
plosively, but this is different from the Big Rip singularity 
in phantom models where the background energy density 
will approach infinity in finite time. Note also that the 
perturbation singularity does not correspond to the cross- 
ing of phantom divide, where l+w^ff changes its sign. In 
some phantom fluid dark energy models crossing this di- 
vide might seem problematical, since fluid perturbation 
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FIG. 6: The effect of the couphng slope on the CMBR and 
matter power spectra. Here fiJii = 0.4 and A — 6.0. Dotted 
hnes are for a = 10, dashed hne for a = 20), and dash-dotted 
for a = 30. The sohd hne is ACDM model 



equations (jB8[) and (|B9[) involve terms proportional to 
1/(1 + w). However, here the possibility of crossing the 
phantom divide has no direct relation to the divergence 
of perturbations linear perturbations. 

As the singularity of the linearized system is ap- 
proached in future in these models, the linear approxi- 
mation certainly breaks down at some point. The per- 
turbative FRW description is no more valid as 5m ^ cxd. 
It seems perhaps plausible that then some kind of matter 
domination is restored, due to the energy creation from 
the Gauss-Bonnet interaction. It might be that some 
singularity really would occur, but this could only be as- 
sessed by considering the field equations beyond the lin- 
ear perturbation theory. Only in the case that nonlinear 
effects would somehow stabilize the growth of ovcrden- 
sities, would the de Sitter phase actually be reached as 
the background calculations imply. Otherwise, it is pos- 
sible that the acceleration is transient. This would be 
favorable from a theoretical point of view, since the S- 
matrix formulation in present versions of string theory 
seems not to be consistent with an eternally accelerating 
universe [t^. 

Previously, instabilities have been found to occur in 



Gauss-Bonnet pre-big bang cosmologv [75ll76j . Such have 
been also recently studied in the context of ghost condi- 
tions in Gauss-Bonnet cosmologies [ttI [TSj . Though as 
well known, when expanded about de Sitter spacetimes, 
ghost modes do not appear in Gauss-Bonnet gravity, it 
has been pointed out that this does not necessarily hold 
in the FRW background since such is characterized by 
non-constant curvature. Quantum field theoretical con- 
sistency requires absence of ghosts, which means that a 
function T multiplying the time derivatives of a variable 
to quantized should be positive. If a function S multi- 
plying the spatial gradient is negative, a related instabil- 
ity might occur when the sound speed Cg = S/T is not 
positive definite. In fact, when an imaginary propaga- 
tion speed appears in an action for a canonical field, one 
would expect already the solutions of its classical evo- 
lution equation to exhibit divergent behaviour. This is 
what we found here. It can be shown that the canonical 
action for the potential $ in the case pm = features an 
effective propagation speed (69| 



-x^[4 + m(5m - 8)] + fi^[6{l + e){n - 1) + y] 

(^_1)[3^2_4(^_1)^2] 



(29) 



When this exceeds one, the speed is superluminal and 
causality could be violated; when it is negative, the sta- 
bility might be lost. Since Eq.([29|l corresponds to prop- 
agation of the scalar field in vacuum, the same Cg would 
be found for any other gauge-invariant variable as there 
is only one scalar degree of freedom. However, in dark 
energy cosmologies featuring the Gauss-Bonnet term one 
should take into account both matter and the corrections 
to Einstein gravity. This does not complicate the anal- 
ysis of the tensor mode, since it decouples from other 
fluids (as long as they are isotropic), and one can just 
plug the background solution (for which matter has been 
taken into account) into an expression for the tensor 
(which is the same as without matter). However, scalar 
modes are non-trivially coupled. In the presence of mat- 
ter, Ea. l|29|) does not describe the evolution of <i>. It might 
still describe the evolution of (50 in some suitable gauge at 
a limit where the impact of matter perturbations on the 
field fluctuations can be neglected, but a priori it is not 
clear that such a limit of cosmological equations exists. 
For the models considered here, one notes by comparing 
Egs. f^^ and pS)) . that the linearized matter perturba- 
tions diverge precisely in the points where the propaga- 
tion speed squared Cg changes its sign. Our interpreta- 
tion is that the fluctuation of the field (p becomes unstable 
when Eq. (|29p dictates, and since it interacts (via gravi- 
tational potentials) with the fluctuation in pm, also the 
latter signals instability at the same instant. Therefore 
the stability of the linear scalar modes is now indeed de- 
termined by a canonical action for any gauge-invariant 
scalar perturbation variable in the vacuum. 
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III. CONSTRAINTS 
A. Cosmological and Astrophysical Constraints 

Armed with all the equations describing the cosmolog- 
ical dynamics, we can now derive the constraints aris- 
ing from astrophysical and cosmological observations. In 
this section we will consider kinematic tests related to 
the background expansion of the Universe: the SNela lu- 
minosity distance-redshift relationship, the CMBR shift 
parameter and the baryon oscillation length scale. 



1. Supemovae la 

The luminosity distance in a flat space is defined as 

dz' 



Dl{z) - (1 + z) 



H{z') 



(30) 



The distance modulus probed by SNela observations is 
then given by 



M 



M = 5 logi 



Dl{z) 
lOpc 



(31) 



We use the "Gold" sample of 157 SNela [3. There the 
observed magnitude m, its error and the redshift is given 
for each supernova. The absolute magnitude M is un- 
known, and its effect is degenerated with Hq entering 
the formula ([31]) from Ea. ([30| . We marginalize over Hq 
(or equivalently, M), by integrating over the likelihood 

as = J x''{Ho)e-^'^"'>'^^^dHo/ J e~^"<-"«^/^dHo. We 
marginalize similarly over the model parameter a. 

For comparison, we perform the likelihood calculatioiis 
also with the Supernova Legacy Survey (SNLS) data [2] 
using then a bit different method. For each of the 115 
supernovae in the SNLS set (labeled here with the index 
i), the distance modulus has been given as a function 
of the stretch factor Si, color factor Ci and the apparent 
magnitude m* as 

7W,;(M, ai, 02) = m* - M + ai{s, - 1) - a2C,. (32) 

We will treat the parameters M, ai and 02 like cosmo- 
logical parameters and find their values which maximize 
the likelihood. For individual supernovae, we keep the 
parameters m*, Si and Ci in their best-fit values. The 
uncertainty in Mi would in principle depend on ai and 
02 , but should be kept fixed while optimizing these global 
parameters 0. Thus we compute 



115 M^{M, ai,a2 
1=1 



5 log 



10 



lOpc 



(33) 



The quantities with lower index ^ we get directly from the 
data, for the parameters M, oi, a2 and /3 we use their 
best-fit values, and for the intrinsic dispersion we use 



cTint = 0.15 for the nearby and aint = 0.12 for the distant 
supernovae. Similar procedure was used in Ref.[7l]. We 
have checked that modifications of the scheme change 
the results only very slightly. The SNLS data set can 
give little bit tighter contours than GOLD set for the 
model at hand. The contours arc shown in dotted lines 
in FIG. [3 

In principle when considering the supernovae observa- 
tions in extended gravity theories, one should take into 
account the possible effects of modified gravity on the 
evolution of the supernovae jsoj. Here we however as- 
sume that the SNela can be treated as standard can- 
dles. If variations in the scalar field (j) can be neglected 
at the scales relevant to supernova evolution, it is de- 
termined by general relativity alone, though corrections 
from Rqb might be crucial at cosmological scales and 
have to be taken into account in the luminosity-distance 
relation (j30p . This assumption needs to be strictly justi- 
fied; we return to a related discussion in Section IIIIB 21 



2. CMBR Peak Locations 

In addition, we compute the CMBR (cosmic microwave 
background radiation) shift parameter ^81,] given by 



n = 



dz 



(34) 



where Zdec is the redshift at decoupling. This number 
TZ captures the correspondence between the angular di- 
ameter distance to last scattering and the relation of the 
angular scale of the acoustic peaks to the physical scale 
of the sound horizon. Its value is expected to be very 
weakly model-dependent, and it can be inferred rather 
accurately from the latest data 84] . We use here the pa- 



rameter values n = 1.70 ± 0.03 and z^ec = 1088+2- The 
resulting contours are plotted with dashed lines in FIG. 



3. Baryon Oscillations 

The imprint of primordial baryon-photon oscillations 
in the matter power spectrum is related to the dimen- 
sionless quantity A in the following way 



hHz,) 



1 

Jo 



dz 



H{z) 



(35) 



The physical length scale associated with the oscillations 
is set by the sound horizon at recombination, which can 
be estimated from the CMBR data ^3]. Measuring the 
apparent size of the oscillations in a galaxy survey al- 
lows one to determine the angular diameter distance at 
the survey redshift. Together with the angular size of 
the CMBR sound horizon, the baryon oscillation size can 
then be a powerful probe of the properties and evolution 
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of the universe. From the large scale structure data one 
can infer that [11] ^ = 0.469 ± 0.017 and = 0.35, as- 
suming the ACDM model. For the scalar spectral index 
we use the best-fit WMAP value 0.95 [sj. The resulting 
constraints are plotted with solid lines in FIG. [9] 



4.. Combined constraints 

The combined constraints arising from the CMBR shift 
parameter, baryon oscillations and the SNela are pre- 
sented in the shaded contours in FIG. [T] There the de- 
pendence on the parameter a is not shown, but we have 
included also likelihood slices on three constant-ri„j slices 
of the parameter space in FIG. ([8]), to explicitly show 
that the value of a is not so tightly constrained, as the 
contours on the a-X planes have are nearly horisontal 
planes. The best-fit values for the are given in Table 
m Note that we have reported the per effective de- 
gree of freedom, Xdo/' ^'^^^ taking into account that the 
coupled model has more free parameters than the ACDM 
model. The fit with CMBR and SNela is as good as with 
the concordance model, but when the constraint from A 
is included the fit gets worse. This conclusion is the same 
for both the Gold and the SNLS data. The latter data 
gives slightly worse fit to the model than to ACDM; when 
7?.-constraint is included the goodness of fit is slightly bet- 
ter than for ACDM; and when the ^-constraint is added 
the goodness of fit is significantly worse, difference of x^ 
being almost Ax^ = 9.5. Taken at face value, the model 
is ruled out by the combining the three kinematical tests. 
The reason why the baryon oscillation data seems to be in 
tension with the other data for this model is that it fixes 
the present matter density rather tightly to the ACDM 
region, ri,„ « 0.27. As a matter of fact in a flat universe 
this constraint is equivalent to requiring that the matter 
density is Vl^ = 0.273 + 0.123(1 -I- wq) ± 0.025, where wq 
is an averaged equation state between the present and zi 
. On the other hand, because the Gauss-Bonnet cou- 
pled field produces a very negative Wef / today, the other 
data sets we included prefer higher matter densities than 
in the standard ACDM model. 

However, on one hand it is not clear how model- 
independent probe of the background expansion the 
baryon oscillation constraint in its present form is. Here 
the the non-negligible contribution of the scalar field to 
the energy density during matter domination probably 
affects the dependence of A on In addition, the 

Gauss-Bonnet modification of gravity might have such 
effects on the nonlinear evolution of galaxy distributions 
that shift the scale at which the oscillation appear us 
today. Finally, in general in dark energy models with 
rapidly varying equations of state the validity of the ap- 
proximation of lumping the observations to the single 
redhift zi has not been establishedp79'|. For such rea- 
sons one has to be careful before ruling out alternative 
cosmologies based on the baryon oscillation constraint. 

On the other hand, as shown in Fig. [51 the baryon os- 
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FIG. 7: The confidence limits corresponding to Ax^ ~ 2.3, 
4.61 and 9.21 for the model in the Q.m - A plane. Dashed 
lines are constraints from the SNela data, solid lines from the 
combined SNela and CMBR shift parameter data, and shaded 
regions include in addition the baryon oscillation scale. Below 
the thick line the scaling solution E is unstable regardless of 
a. In the upper panel the SNela data is the Gold data set, 
and a (together with Ho) is integrated over 1.5A < a < lOA. 
In the lower panel the SNela data is the SNLS data set and 
a (together with M, a\ and 02) is set to minimize the . 



dilations and the CMBR acoustic scale alone or together 
do not limit the allowed parameter space very strictly. 
The best fit to these observations is x^ ^ 1; and they 
are compatible with lower matter densities than with the 
SNela data included. As mentioned previously, it is pos- 
sible that the supernova constraints do not apply in their 
present form to the model, and hence it is interesting to 
note that when they are discarded, the model matches 
excellently with the other observations. It is then also 
consistent with high values of A, and the tension with 
nucleosynthesis bound, to be discussed next, disappears. 



5. Nucleosynthesis 

The nucleosynthesis bound restricts the allowed 
amount of dark energy in the early universe. If the ex- 
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FIG. 8; Parameter constraints in the {a/\, A) plane. In the 
right panels, dashed contours (corresponding Ax^ ~ 2.3) are 
constraints arising from SNela and the solid contours include 
also TZ. In the left panels, we have in addition included the 
A. There at the solid lines one has G', = at the present. 
In the upper pair of figures f2m = 0.3, in the middle panels 
= 0.35 and in the lowest two figures Q.m = 0.4. One 
notes that the background expansion does not very sensitively 
probe Q, and that requiring = restricts very strictly the 
parameter combinations. 
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TABLE L The best-fit values for ACDM model compared 
with fits of the coupled scalar field model for some parameter 
values when a is marginalized over. The error limits corre- 
spond to change of unity in the . The degrees of freedom 
in the first row are 157 — rf, in the second 158 — d and in the 
third 159 - d, where d = 1 for the ACDM and d = 2 for the 
G0 models. 



pansion rate is not due to radiation alone, the prediction 
for the abundance of the light elements produced during 
nucleosynthesis will be modified. This places constraints 
on Q.^ &i a ^ 10~^°. As tight constraints as fi^ < 0.045 
have been reported [13], but more a conservative limit 
is < 0.2 0. Here one should keep in mind that the 
observationally preferred value oiujbh? might turn out to 
be different than in the concordance models (as it seems 
to be with ajmh"^- Nevertheless, if the scalar field is track- 
ing during nucleosynthesis, the more conservative bound 
translates into A > 6.3. While with potential slopes cor- 
responding to such values, a good fit can still be achieved, 
the best fit values are typically at smaller A if the SNela 
data is taken into account. Therefore one might want to 
consider relaxing the nucleosynthesis constraint. 

There are several ways to do that. The simplest op- 
tion is to assume that the scalar field has not reached the 
tracker solution yet, but begins to approach it during or 
after the nucleosynthesis. Then one just considers large 
enough initial values for the field, so that the potential 
term in Eq.dH]) is small compared to the Hubble drag at 
very early times. This could be the natural outcome of 
inflation [85j. If one however insists that the field must 
have entered the tracking phase well before nucleosynthe- 
sis, it would be necessary to modify the early dynamics 
of the field. If the slope of the potential is steeper for 
small values of the field, the energy density residing in 
the quintessence field is subdominant until it reaches the 
shallower region, and choosing the parameters suitably, 
this could happen after a sufficient amount of light ele- 
ments is produced. It might be reasonable to consider A 
as a function of the scale factor [sst. Alternatively, an ad- 
ditional field dependent features in the coupling / could 
be used to hold the field subdominant enough in very 
early universe. 



FIG. 9: The confidence limits corresponding to A^^ = 2.3, 
4.61 and 9.21 for the model in the Qm - A plane. Dashed lines 
are constraints from the CMBR shift data alone, solid lines 
from the baryon oscillation scale alone, and shaded regions 
combine these data. The coupling a (marginalized over, in 
the range 1.5A < a < lOA. 



B. Solar System Constraints 

Time variation of the gravitational constant is tightly 
constrained by observations. These observations include 
various tests of the force of gravity within the Solar sys- 
tem and laboratories, and indicate that {dG^,/dt)/G^ is 
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less than about 10 per year, where G* is the effective 
gravitational constant 58]. This bounds translates into 

|g^|<0.01. (36) 

To derive the variation of this constant for the coupled 
Gauss-Bonnet gravity, we follow the approach of Ref. [s^ 
where cosmological perturbation equations were consid- 
ered at their Newtonian limit, i.e. assuming small scales 
and small velocities the Poisson equation was derived, 
and the effective strength of gravitational coupling was 
read from the resulting expression, which relates the gra- 
dient of the gravitational potential to the perturbations 
in matter density. We will first compute results ensuing 
from this approach and then discuss the Post-Newtonian 
parameterization. 

1. Newtonian limit 



5 




a 



FIG. 10: Variation of the gravitational constant, G^/G, at 
the present as a function of the coupling a. The solid line is 
for the first model in Table HIl and the dashed for the second 
one. 



Consider subhorizon scales at a limit where the time 
derivates of the perturbations are much smaller than 
their gradients. In the Newtonian gauge we have ip = 
and ipA = ^- Then k = 3H'i>. This gauge is also called 
zero-shear or longitudinal for the reason that there x = 0- 
The ADM energy constraint Eq. (|B2p is then, at this limit 

^ = -2(1 - 8Hf )<P - SH^Sf. (37) 

The shear equation (|B4p reveals there is effective 
anisotropic stress in the sense that $ 7^ — 'F: 

(1 - 8i?/)* + (1 - 8/)$ = -8{H + H^)6f. (38) 

To obtain the effective Poisson equation here, we need 
to eliminate the scalar field perturbation from the above 
two equations. The fluctuation Scj) is given by the Klein- 
Gordon equation (jB6p . 

Scj} = 8f[2{H + H^)<^> - H'^^]. (39) 

From the previous three equations, it is then straightfor- 
ward to obtain the Poisson equation, 

^3-* = 8nG,pS. (40) 

We find that the effective gravitational constant here co- 
incides with the expression (j25p . This is in spite that 
the Newtonian limit is considered at a static configura- 
tion of the gravitational sources. This takes the approx- 
imation we used in Section Bp further, since there we 
derived the small scale limit of the dynamical equations. 
We can undertand that the two G* agree since the time 
derivatives of say ^E" are expected to be proportional to 
\]/ H'if. The gradient terms are more important at 
small scales and hence determine there the main contri- 
bution to the gravitational coupling of evolving struc- 
tures. 



Though not obvious from the formula (PS)) . it equals 
one when the coupling goes to zero. Generally, when 
the coupling is significant, i.e. fi is of order one, then 
one expects the G'^/G^ to be of roughly of order one as 
well. This is confirmed by numerical evaluation of the 
time variation of G* for the model considered here"^. As 
claimed in Ref. [86j, one has to assume "an acciden- 
tal cancellation" to satisfy the bound Eq. ([M]) in the 
presence of significant Gauss-Bonnet contribution to the 
energy density. This means that we have to fine-tune 
the coupling a in order to eliminate the time variation 
of the effective gravitational coupling. Then the Newto- 
nian limit in general exhibits time- varying G* , but at the 
present this G* appears to us to be a constant. In prac- 
tice, the requirement of well-behaved Newtonian limit 
sets the magnitude of the coupling a. Such a limit does 
not exist for all A and flm- When it is found and it re- 
quires a to be fixed with the accuracy of about 0.01, see 
FIGs. [5] for actual constraints and [TO] for the variation as 
fimction of a. 

In Table |IT] we report the Xdof ' computed under the 
condition G* = 0, for some parameter choices for our 
model. Even with the coupling fixed in to yield G* = 0, 
the coupled Gauss-Bonnet quintessence fits both the 
SNela and the CMBR shift parameter data as well as, 
or slightly better than the ACDM concordance model. 
However, again this is only when the baryon oscillation 
constraint is not taken into account: with the A com- 
bined with all the other data sets including the Solar 
system constraints the model is in blatant contradiction 
with the data. 



^ It is possible to find an expression for the G'^/Gt in terms of x, 
y and fi, but it is somewhat lenghty. 
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Data set 


ACDM 


R^g<f) model 




Xdof 




Xdo/ 




A 


a 


SNela 


1.142 


0.314 


1.146 


0.42 


5.1 


32.3 


SNeIa+7^ 


1.144 


0.277 


1.141 


0.44 


5.2 


33.8 



TABLE II: The best-fit values for ACDM model compared 
with fits of the coupled scalar field model for some parameter 
values when the coupling a is set in order to fix the present 
time variation of G* to zero. The degrees of freedom in the 
first row are 157 — d, in the second 158 — d and in the third 



159 — d, where d - 
models. 



1 for the ACDM and d = 3 for the Rcb(, 



inadequate and setting a to a certain value would not 
guarantee the viability of the model). Such an adjust- 
ment would apparently ruin the potential of the model 
to alleviate coincidence problem, since the slope of the 
coupling function f{(j)) would change just today. How- 
ever, it might be possible to associate the change of the 
slope with the triggering of the acceleration, since it is 
then that the curvature corrections become dynamically 
important and one could expect higher order modifica- 
tions to enter in the play. We hope to address these issues 
more quantitatively elsewhere. 



2. Schwarzchild metric 

The rather tight bounds we get could be loosened if 
one takes into account that cosmological variations of G 
and other gauge-couplings might be different from the 
ones we measure on Earth or within our Solar system 
[stI . [ssI . [sol . [90| . The model is set up in such a way that 
the corrections to the Einstein gravity will affect the over- 
all expansion of the Universe, while the Solar system is 
clearly subcosmological in scale and nonlinear in nature. 
Thus the local values of the scalar field and the potentials 
could be something else than at cosmological scales where 
the linear perturbation equations of Appendix Eq.([B]) are 
expected to apply. If for example, the scalar field hap- 
pened to be nearly frozen at our neighbourhood we would 
observe G* = G though at larger scales the linear pertur- 
bation would evolve according to Eq. where possibly 
G', ^ 0. 

Conventionally the constraints ensuing from Solar sys- 
tem experiments are imposed on the deviations from 
the Schwarzchild solution written in form of the Post- 
Newtonian parameters (PPN). Here we computed the 
Newtonian limit of the cosmological perturbation equa- 
tions in the FRW background to derive the constraints, 
and there is reason to doubt that results ensuing within 
the PPN formalism would not be equivalent. Therefore 
our considerations of the Newtonian limit in the previ- 
ous subsection must be regarded as preliminary. More 
detailed study of the Solar system constraints is left for 
for future. Since the is quadratic, its value for the 
Schwarzchild metric goes like However, due to the 

dynamics of the field, this does not necessarily mean that 
the curvature correction would dominate at small scales. 
Previously local constraints for a scalar-Gauss Bonnet 
coupling have been investigated in the case that the po- 
tential can be neglected 91, 92]. Even then all the effects 
of the scalar field can be negligible for any Solar sys- 
tem experiment provided f"{(t>o) is negative. One also 
notes that due to specific properties of Rq^ , if the cou- 
pling f{4>o) happens to be in its minimum at present all 
the corrections trivially disappear. With these remarks 
in mind, it seems easy to adjust the function f{4>) in 
such a way that the local constraints would be satisfied 
(even if our preliminary analysis would turn out to be 



IV. CONCLUSIONS 

We studied cosmological phenomenology of dark en- 
ergy based upon a low-energy effective string theory ac- 
tion featuring a compactification modulus and taking 
into account the leading order curvature corrections. One 
is then lead to consider a generalized scalar-tensor theory 
which includes a coupling to the Gauss-Bonnet invariant. 
Specifically, we investigated the evolution of perturba- 
tions. There the simple closed form equation for the 
linear matter inhomogeneities, together with its effective 
gravitional constant ((25)) could be highlighted as one of 
the main results of the paper, since it enabled to find sev- 
eral new model-independent results of perturbation evo- 
lution in Gauss-Bonnet dark energy, these including the 
following. 

• The evolution of matter perturbations is scale- 
invariant at small scales in the presence of the 
Gauss-Bonnet term, and thus the shape of the mat- 
ter power spectrum is retained. 

• The growth rate of matter perturbations, which is 
easily extracted from Eq.p4|). can be compatible 
with observations even in the presence of significant 
contribution from the Gauss-Bonnet interaction. 

• The stability of perturbations can be read off from 
the expression for the effective gravitational con- 
stant. A divergence is possible. 

This gravitational constant cannot be deduced by match- 
ing Schwarzchild and FRW -type metrics for subhori- 
zon evolution of spherical matter overdensities, as the 
Jebsen-Birkhoff theorem is not in general respected in 
these models. The Newtonian limit, as discussed in Sec- 
tion |IIIBl does feature this effective strength of gravity, 
which might be used to efficiently constrain the coupling. 
However, the relevance of this limit to the experiments 
within Solar system is an open question to be addressed 
more carefully in forthcoming studies. 

In addition to such general considerations mentioned 
above, we chose a particular model and scrutinized its 
predictions for cosmology. Specifically, we characterized 
the potential and the coupling by exponential functions, 
as with the Nojiri-Odintsov-Sasaki modulus This 
parameterization is appealing in its minimality, as within 
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it one can describe Gauss-Bonnet dark energy with just 
two extra degrees of freedom compared to the ACDM 
model. Moreover, the exponential forms are well moti- 
vated on fundamental grounds, and allow dark energy 
solutions without introducing unnatural scales into the 
Lagrangian. Many of the previous cosmological studies of 
this type of low-energy string theory have been concerned 
with asymptotic solutions and confined to cases where 
the Gauss-Bonnet term is subdominant, and thus only in- 
troducing small deviations to standard quintessence sce- 
narios. Here allowed a crucial role for the Gauss-Bonnet 
term and examined its consequences in detail. Such ap- 
proach also proved to be very useful, as it revealed several 
possibilities of cosmological evolution present already in 
this simple model, these including the following. 

• The Gauss-Bonnet coupling can act to switch the 
decelerating expansion into acceleration after a 
scaling matter era, thus perhaps alleviating the co- 
incidence problem. 

• The curvature interaction may momentarily push 
the universe into a phantom era, which will then 
not lead into a Big Rip singularity in the future. 
Hence the model can explain superacceleration. 

• However, the linear perturbations might diverge 
in the future. This might terminate the dc Sitter 
phase, possibly helping to make contact with string 
theory. 

Let us briefly comment each of these possibilities. 

The field is drawn into the scaling attractor virtually 
regardless of initial conditions, and at late times accel- 
eration can be onset provided simply a > A. This can 
be seen as a possible way towards solution of the coinci- 
dence problem, though it is occasionally considered that 
in a preferable solution the future evolution of the Uni- 
verse would be characterized by a constant ratio of fi^ 
and r2„i different from 1. However, a transition from a 
scaling matter era to an accelerating scaling era has been 
shown to be impossible for a very general class of scalar 
field models [9^. Another remark is that the question 
of why dark energy is beginning its domination today 
is now translated into the question that why the cou- 
pling comes into play just recently. The same holds for 
any other coupling mechanism proposed so far, but it is 
still fair to say that such mechanism can provide an ap- 
proach to tackle the coincidence problem. The relevant 
dynamics are not always obvious when only the asymp- 
totic behaviour is considered, and therefore it was only 
very recently j3§, |41J found that the simple model con- 
sidered here can (without flipping the fleld to a phantom, 
i.e. putting 7 = — 1 in the action ^) feature a transition 
from a matter dominated to an accelerating phase. 

We also showed that during this transition, the dy- 
namics of the coupling can cause a phantom expan- 
sion. Again, this phenomenon is unaccessible when only 
asymptotic behaviour or vacuum is taken into account. 
Several authors have shown that phantom dark energy 



does not necessarily lead a Big Rip singularity in the pres- 
ence of the Gauss-Bonnet term. Here we found that this 
term can in fact ally with a canonical scalar field into an 
effective phantom energy. This is clearly a more appeal- 
ing way to produce Weff < — 1, since thereby the phan- 
tom era can take place without the introduction strong 
energy condition violating components, while the Big Rip 
singularity is still avoided. 

However, one of the peculiarities of this class of models 
is that their linear perturbations can diverge, even while 
the background solution is not singular. We argued that 
such a possibility does not necessarily indicate an incon- 
sistency of the model. For the specific model considered 
here, the divergence typically awaits in the future when 
the evolution would begin to asymptote to the de Sitter 
expansion. The dynamics is such that the (linearized) 
inhomogeneities grow explosively as the ssc < limit 
is approached. Hence the linear approximation perhaps 
together with the whole FRW description breaks down. 
It is tempting to believe that the de Sitter phase will 
not then be reached but the inhomogeneities take over. 
This would enable to reconcile the accelerating universe 
with string theory in these models in yet another way. 
Since if the acceleration is transient one can consistently 
define the set of observables in string theories [Tij. Their 
present S-matrix formulations seem to be in odds with an 
eternally accelerating universe, like in the concordance 
model with the cosmological constant. Would the uni- 
verse however somehow pass through c| = 0, one would 
have to invoke higher-order curvature corrections or mod- 
ifications of the exponential parameterization to escape 
the possible singularity. 

We calculated the implications of Gauss-Bonnet dark 
energy to various different cosmological and astrophysical 
observations and matched them with several data sets. 
The general conclusion we arrive at is that the quadratic 
curvature coupling has interesting and non-trivial sig- 
natures while being in a good agreement with a large 
amount of data. We checked how the Gauss-Bonnet in- 
teraction affects the evolution of scalar field fluctuations 
and how these in turn impact the large scale structure 
and CMBR observations. Using combined datasets re- 
lated to the CMBR shift parameter, baryon oscillation 
scale, supernovae la luminosity distance and variations 
of the effective gravitational constant in the Solar sys- 
tem, we derived confidence limits on the parameters of 
the model. The study presented not only shows the re- 
quirements to be imposed on models of Gauss-Bonnet 
dark energy for them to exhibit completely realistic cos- 
mology, but also quantifies the degree to which these re- 
quirements can be within the most minimalistic set-ups. 
Due to its simplicity, matching with observations, capa- 
bility to shed light on the coincidence problem, possi- 
bility incorporate phantom (and perhaps also transient) 
acceleration, and not least its theoretical motivation, the 
Gauss-Bonnet quintessence seems to be a promising al- 
ternative to a cosmological constant. 
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We end by mentioning some issues yet to be pursued. 
Among these are the nature of the divergence and the 
Solar System hmit in these models. Also, we found 
that there is some contrast with the amount of early 
quintessence and the nucleosynthesis bound and with the 
present matter density and the baryon oscillation scale. 
Whether these enforce one to resort to more elaborate 
models depends on several assumptions, like the validity 
of the baryon oscillation constraint and the initial condi- 
tions of the scalar field after inflation. One caveat is also 
that these contrasts can disappear if the Gauss-Bonnet 
modification affects the intrinsic evolution of supernovae 
and thus changes the effective luminosity-distance rela- 
tionships of these objects. This does render some of our 
conclusions indecisive, but in an interesting way. Sup- 
pose that the coupling Gauss-Bonnet modification does 
influence small scale physics considerably. Then one has 
to consider f{(j>) such that it evolves to a constant as 
dark energy domination begins, so as to comply with 
the Solar system experiments (which is very much the 
approach in standard scalar-tensor theories). Then one 
could match all the observations, while the expansion 
of the universe at small z could be something consider- 
ably different than commonly established from the SNela 
luminosity-redshift curves. 
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Here w is given hy w = n°/{3{n°+an'i,)). Since the field 
is attracted to the scaling solution E regardless of the 
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of the coupling and the potential, the initial condition for 
jj, is determined by the required amount of matter today 
The fractional matter density flm as a function of 
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APPENDIX A: BACKGROUND 

In the appendices we will consider some technical de- 
tails relevant in solving the cosmological equations. The 
first appendix concerns only the kinematic evolution 
model specified by Eq.Q. The equations in the second 
appendix are for perturbations in general for an action 
In fact, when one sets 7 = 0, they describe also 
modified Gauss-Bonnet models [i^ featuring a function 
fi^Gs) their action, since there exists a simple map- 
ping between the theories p3|. 

Employing the dimensionless variables (|10p we can 
write the system (I13|) to derive the explicit evolution 
equations for each individual variable. It would be pos- 
sible to algebraically solve all but two of them, but we 
rather use the three variables related to the fractional 
energy densities of the kinetic term of the scalar field x, 
its potential y and the coupling fi. The coupled ordinary 
differential equation system for these three variables is 
then 



(Al) 



(A2) 



(A3) 



(A4) 



log(a) is given by f7,„ = (3 — — 3/i- — y)/3. Since on the 
other hand — ^^a^^^^^"^^ / , we can readily ob- 
tain the Hubble parameter as function of redshift from 
the solution for x, y and fi. One consistency check that we 
have made is the comparison of the numerical derivative 
of the Hubble parameter with the algebraic solution for e, 
Eq. (|A4|1 . One might be concerned that the system (|Aip - 



, _ 2e(-3 + x^ + 6fi) + 3.x^(l + w) + 3(1 + w)(-3 + 3n + y)+ x{3afi - \y) 
^ ^ 2a;2 + 3fJ. ^' 

y' = -{2e + xX)y, 

, _ 2ax^ + 2e(3 + x^ - 3fi) - 3x^(1 + w) + xXy - 3(1 + w){-3 + 3/i + y) 
^ ~ 2x2 + 3/x 

where the slow-roll parameter e is given by 

2ax^H - 3^*2 + 2a;''(-l + w) + 2x'^[^j.{-1 + 3w) + (1 + w){-3 + y)] + Xx^J.y 
^' " 4a;2(-l + - 3^2 ■ 

I 
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is not well behaved in the limit that the denomina- 
tors in the equation go to zero. The denominator vanishes 
when the field is constant. However, then also the numer- 
ators go to zero. We have found no numerical problems 
in the situations when the sign of the field velocity hap- 
pens to get flipped. On the other hand, the denominator 



of Eq. (jA4p would vanish also when ^ 



2^2 



This would occur only when the field would be decreas- 
ing at enormous speed, and such a configuration is indeed 
impossible to reach. 

For completeness, we also mention that for the fixed 
point G we have 



10 



100 - 



1 



1000 - 54a2 + 9ri/3 3a 



(A5) 



and 



A* = 



27\'^ + irk + lOO/c^ + 20(500 + 3r + SO/c) 



3a2(320-f 26fc-|-5fc2) /(9fc 



(A6) 



where r = a(2000 - 2Ma^ + Qa^ fl'^ and k ^ (1000 



54a2 -f 9r)i/3. 



APPENDIX B: LINEARLY PERTURBED 
EQUATIONS 

In this appendix we consider the scalar perturbation 
equations for coupled Gauss-Bonnet gravity in a general 
gauge. These equations are consistent with the ones pre- 
sented in Refs.Kmlsgl. 



The metric is defined as in Eq. (fT6|) . 
variables satisfy the constraints 



Recall that the 



X = al3 + a^^, K = 3{HipA-'p) + ^X- (Bl) 



Due to general structure of covariantly modified gravity 
[66l |. the last three conservation equations are deducible 
from the set of field equations (IB2|) - (|B5|) when the matter 
content is known. 



The energy constraint (Gq component of the field equation) is 



(l - 8Hf) 



7^- 



1 



5p,r, + i4>H + V'54>) + 4i/2 iH6f + —5f 02) 



and the momentum constraint {Gf component) is 



l-8i// 



fc2 



a{Pm +Pm)Vm + 



12H^ 



(B3) 



The shear propagation equation (G* — I'^JC^ component) reads 



(l - 8i//) X 



H-8{Hf + Hf + H^f 



X - (l - 8iJ/) ^A-(l- 8/) ifi - n™ - 8 (i/ + H^) Sf. (B4) 



The Raychaudhuri equation (G^ — Gq component) is now given by 



l-8Hf]k + (2H -8Hf -8Hf -12H^f] K-12H^fipA 



1.2 

3H-^ + 2702 _ 24H^f - A8HHf - UH^f + 8Hf— ] ipA 



i((5p,„-K3(5p,„) + 270,5'(^~y'(5 



SH'^Sf + (2H + H'^^ (^3HSf + ^Sf 



(B5) 



The previous three equations can be considered as evolution equations for ip, x ^-^id '^i respectively. For instance, the 
Raychaudhuri equation might be written as k = sources/(l — 8Hf). Also the other equations imply a divergence of 
perturbations in the case that 8Hf = 1. In the models considered here, 8Hf < 1 is always satisfied; this condition in 
fact is necessary for absence of tensor ghost modes [zl]. Similar considerations apply for the factor (1 — 8/) appearing 
in some of the equations. One might understand this by noting that the first condition for stability of scalar modes is 
automatically satisfied given the conditions for stability of the tensor modes [78j]. The second condition for stability 
of scalar modes, which is not always guaranteed in models considered here, is discussed extensively in Section [Tl B 31 
The equation of motion for the scalar field is 



■ + SHd4>) + ( — + y" ) 5,/) = 7</) {ipA + k) + 7 



(24> + 3H,p 



^A~I"RlBH-f5RlB, 



(B6) 
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where 



SRl;j^ = 4:H^6R - 16H {hk - , 5R^2 

The continuity equations for matter are 

5m = 3H (Wm - C^) (^m + (1 + Wm) 



(B7) 



(B8) 



1 Hv„ 



1+Wm 3{1 + Wm)a'^Pr, 



(B9) 



To implement the model in the CAMB code, we have included there an integration of the background equations 
described in the Appendix [X] The solution is then interpolated to obtain it at arbitrary points. The perturbation 
equations are modified to be consistent with the set of equations listed above, as written in the synchronous gauge 
(with the CAMB variables slightly different from those h, rj and x discussed in III Bp . Modifications have to be 
taken into account in the evolution equations for perturbations to obtain the matter power spectrum, and in the 
evaluation of their line-of-sight integral to obtain in additio n the CMBR spectrum. 
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